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We show that a thermoelectric transport through a Quantum Dot (QD) - single-mode Quan¬ 
tum Point Contact (QPC) nano-device demonstrating pronounced fingerprints of Non-Fermi Liquid 
(NFL) behaviour in the absence of external magnetic field is protected from magnetic field NFL 
destruction by strong spin-orbit interaction (SOI). The mechanism of protection is associated with 
appearance of additional scattering processes due to lack of spin conservation in the presence of both 
SOI and small Zeeman field. The interplay between in-plane magnetic field B and SOI is controlled 
by the angle between B and Bsoi • We predict strong dependence of the thermoelectric coeffi¬ 
cients on the orientation of the magnetic field and discuss a window of parameters for experimental 
observation of NFL effects. 

PACS numbers: 73.23.Hk, 73.50.Lw, 72.15.Qm, 73.21.La 


I. INTRODUCTION 

The paradigm of Landau Fermi Liquid (FL)^ is one 
of the cornerstones of modern condensed matter theory. 
Based on the concepts of quasiparticles - well defined ex¬ 
citations whose energy in the long-wave limit is greater 
than their decay rate, the FL theory successfully explains 
the behaviour of normal and superconducting metals giv¬ 
ing universal predictions for thermodynamic and trans¬ 
port propertied. The FL phenomenology is justified in 
many microscopic models describing interacting fermions 
in- and out- of equilibrium. However, there are several 
cases where a violation of the FL picture is observed ex¬ 
perimentally (e.g. in strongly correlated electron sys¬ 
tems such as heavy fermion compoundd, unconventional 
superconductors 2 and quantum transport through nano¬ 
structures 4 5 ). The pronounced Non-Fermi Liquid (NFL) 
behaviour of these systems is attributed to a breakdown 
the quasiparticle concept: the decay rate of low-energy 
excitations becomes greater than the energy of the exci¬ 
tations itself. 



FIG. 1. (Color online) Typical setup for thermoelectric mea¬ 
surements: ’’cold” contact (light orange area) at reference 
temperature T - Quantum Dot - Quantum Point Contact 
electrostatically defined by gates (blue boxes), separated by 
a tunnel barrier from a ’’hot” (deep orange) contact at tem¬ 
perature T + AT. The voltage V is applied across the device 
for zero current measurements (see text for the details). 


While the fingerprints of NFL physics in thermodynam¬ 
ics of strongly correlated systems and quantum transport 
had been seen experimentally^, it is also generally ac¬ 
cepted that the NFL picture is extremely sensitive to 
variation of external parameters being unstable against 
the FL ground state. Thus, the stability of the NFL 
domain and the possibility to observe strong deviations 
from the Landau FL paradigm posses major challenges 
including the development of theoretical models and pre¬ 
dictions for the stabilization of NFL-states. Important 
questions are: is it possible at all to protect unstable 
NFLs? What are the physical observables which demon¬ 
strate the most pronounced manifestation of the NFL 
physics? 

In this paper we present an example based on a win¬ 
dow of parameters within which the observation of strong 
deviation from the FL picture can be protected and ex¬ 
tended by effects of strong spin-orbit interaction (SOI). 
The physical observables we consider in this work are 
thermoelectric coefficients of a nano-device (Fig. 1). Our 
theoretical model justifying NFL behaviour is a two- 
channel Kondo (2CK) model 6 ^. While the scattering 
of single orbital channel electrons on a resonance quan¬ 
tum impurity itself leads to strong modification of the 
thermoelectric transport properties within the Landau 
FL paradigm through strong renormalization of the FL 
energy scalded], the detour from the FL picture is pre¬ 
dicted to change co mplete ly both electricR 2 15 and ther¬ 
moelectric transports 16 17 . For example, one of the man¬ 
ifestations of the NFL behaviour in quantum transport 
is associated with the logarithmic enhancement of the 
thermoelectric power 16 in the situation when the 2CK 
model originates from the charge Kondo effect in a sin¬ 
gle mode quantum point contact (QPC) - quantum dot 
(QD) setup tuned by gate voltages to the Coulomb block¬ 
ade (CB) peak regim^M^. j n that case two channels 
are the electron spin degrees of freedom while the al¬ 
most transparent QPC (weak back-scatterer) works as a 
quantum impurity. The 2CK physics is known to be un¬ 
stable with respect to any effects which can potentially 
break (statically or dynamically) the symmetry between 
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the channels 5 17 . In particular, it has been shown that 
the effects associated with time-reversal symmetry break¬ 
ing (TRS) due to an external magnetic field restore the 
FL properties at temperatures below T e ff tunable by the 
field 17 . The universality class of unstable 2CK model 
than changes to single channel Kondo problem (1CK). 
Fully screened 1CK is characterized by stable local FL 
properties. Therefore, while being very attractive from 
theoretical point of view, the 2CK physics suffers from 
serious experimental obstacle^ 4 5 18 impeding a direct ob¬ 
servation of NFL behaviour. 

The paper is organized as follows: we describe possi¬ 
ble experimental setup for observing the NFL transport 
in Section II. A theoretical model accounting for the in¬ 
terplay between SOI, external magnetic field and effects 
of Coulomb blockade in Quantum Dot is presented in 
Section III. We discuss the solution of one-dimensional 
quantum-mechanical scattering problem in the presence 
of strong SOI in Section IV. An effective model describ¬ 
ing a low-energy physics of the problem and its exact 
solution is presented in Section V. The transport coeffi¬ 
cients computed with the help of the exact solutions are 
discussed in Section VI. Section VII is devoted to discus¬ 
sion of the key results of the paper including estimation 
for the parameters and definition of the conditions nec¬ 
essary for experimental observations of the NFL physics. 
Summary and Conclusions are given in the Section VIII. 
Details of derivation of the effective model are presented 
in the Appendix A. 


II. PROPOSED EXPERIMENTAL SETUP 

We consider a two-terminal nano-device (see 
Fig. 1) d esign ed to be used for thermoelectric 
measurement f®—. The QD - QPC contains 2-d 
electron gas (2DEG) confined in ^-direction (light 
orange area on the Fig. 1). The open QPC connects it 
to the drain at the reference temperature T. We assume 
that the Rashba SO l 21 * 22 1 (caused by the gradient of 
the confining potential in ^-direction) leads to appre¬ 
ciable effects which we will discuss in this paper. The 
source is separated from the QD by a tunnel barrier 
with low transparency |t|<Cl. The temperature of the 
source (deep orange) is adjusted by the Joule heat 
controlled by the current Ij flowing along the lead 
(black arrow). The temperature difference AT across 
the tunnel barrier is assumed to be small compared 
to the reference temperature T to guarantee the linear 
response operation regime for the device. The QD is 
electrostatically controlled by two plunger gates (blue 
rectangles) to adjust the size of the electron island. The 
device is operated in the steady state of zero source-drain 
current I sd =0=GAVth-\-GTAT, controlled by applying a 
thermo-voltage A V t h between the source and the drain. 
The QPC (denoted by the cross in the light orange area) 
is tuned to the single mode regime characterized by 
a controllable small reflectivity |r|<Cl. Under this as¬ 


sumption and neglecting the resistance of the ” metallic” 
QD we assume that the voltage difference A Vth arises 
across the tunnel barrier between the source and QD. 
The transport coefficients: electric conductance G and 
thermoelectric coefficient Gt (measured independently) 
define the thermoelectric power (TP) S: 


Gt 


dhd r = 9I sd 
<9 A T’ 8V 


AV th 


AT 


I s d = 0 


Gt 

~G m 


We assume that the magnetic field (blue arrow) is applied 
parallel to the plane of 2DEG to avoid orbital effects. 


III. THEORETICAL MODEL 


The theoretical description of setup (see Fig 1) is for¬ 
mulated in terms of the Hamiltonian: 


H = H S + H d + H tun + H z . (1) 

Here H s and H d are the Hamiltonians of the source 
(’’hot” contact) and the drain (’’cold” contact), respec¬ 
tively. H tun describes tunneling between the source and 
the drain and H z accounts for the Zeeman effect in both 
contacts. We assume that the source can be described 
by a standard FL approach 

Hs = ^ ^ e ka c \ a c ka ( 2 ) 

k,a 

here c t and c are creation/annihilation operators of quasi¬ 
particles (we adopt a system of units h = kp = 1). The 
drain H d = H c + Hqpc includes the Coulomb blockaded 
QD described by charging Hamiltonian H c and QPC rep¬ 
resented by 

Hqpc = Hq + Hsoi + Hps (3) 


We assume that the charge Q = e(n s + n d ) in the QD 
is weakly quantized (mesoscopic CB regime 27 ) and con¬ 
trolled by the gate voltage V g : 

H c = E c [h s +n d -N(V g )] 2 , (4) 


here n s and n d are the operators of the number of elec¬ 
trons that entered the dot through the source and the 
drain respectively, E c ~ e 2 /Lqd is the charging energy 
of QD with geometric size Lqd. Below we ignore ef¬ 
fects associated with finite mean-level spacing in the dot. 
While charge is only weakly quantized in the mesoscopic 
CB regime, the spin remains a good quantum number in 
the absence of SOL However, when the SOI is present, 
two spin sub-bands are split horizontally in k-space and 
while spin is no more conserved, the sub-band index 
characterizes quantized states instead. The single mode 
QPC being a short quantum wire can be viewed as 1- 
d electron system in the presence of Rashba SOI 23 26 


Hsoi = 


k x n z \ -cr: 


}■' 


H 0 = —ivp 



(: y ) dy^x,a (: y ), 


( 5 ) 
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Hsoi = olr^f ^ A 

Act 


/ °° 

dy + ( 6 ) 

-OO L J 


We denote here by ^a,ct the left (A=—) and right (A=+) 
movers with spin a =t,l- The constant qlr charac¬ 
terizes Rashba SOI strength. The kp and VF=kF/m* 
correspond to the Fermi momentum and Fermi velocity 
(here m* is a fermion’s mass). The 1-d electron trans¬ 
port through the QPC is along the y -axis (see Fig.l). 
The Rashba SOI Hsoi=&Rk y cr x is associated with the 
electric field gradient along the z -axis and can be char¬ 
acterized by the effective SOI field gyB^soi/^—^ r^f^x 
perpendicular to the direction of electron transport (g is 
the Lande factor, fip is the Bohr magneton). Notice, that 
the SOI field alone does not lead to the TRS breaking. 

The backscattering (BS) Hamiltonian describes a scat¬ 
tering of electrons with momentum transfer 2k f on a 
non-magnetic quantum impurity located at the origin 
and characterized by a short-range potential V(y ): 

H B s = Y,j dy^l a (y)V(y)^- Ka (y)e- 2iXkFy . (7) 

A ,cr 


The Hamiltonian H tun represents the weak tunneling 
I tk | = \t \"Cl of the electrons from the left contact to QD: 


H tun 


kXa 


[tkcl^Xai-Oo) + h.C. 


(8) 


The Zeeman Hamiltonian H z describes the effects of the 
external magnetic field H z = —g/iBE{s s + s^), where 
s s and Sd are the spin densities of electrons in the 
source and drain respectively. We consider a situation 
when both sizes of the QD ( Lqf> ) and QPC ( Lqpc ) 
are small compared to the SOI length scale Lqf> ~ 
LQpc^lsoi = l/(m*&R)- Since the effective energy scale 
determining the behaviour of the transport coefficients of 
the model (| 2 ][ 8 | which will be referred below as the Kondo 
temperature T K isP ^ E c (see Appendix A), the con¬ 
dition Isoi^Lqd is equivalent to gfiBBsoi^TK (see a 
discussion about interplay between Kondo effect and SOI 
in Ref. 28). We also assume that the SOI effects in the 
QD are already taken into account by using the approach 
developed in Ref. 29). 


IV. SCATTERING IN THE PRESENCE OF SOI 
AND MAGNETIC FIELD 

We consider Id scattering problem in the presence of 
SOP and Zeeman field applied parallel to the plane of 
2DEG. The Hamiltonian is given by 

U = 'Ho + V{y) = -T- + a R a x k -7 a-B + V(y). (9) 
2m* 

The short-range potential V( y) des cribes a non-magnetic 
impurity located at the origirPI 2 S[ The electron’s trans¬ 
port is along ^/-direction, k=k y . Angle cp characterizes 



FIG. 2. (Color online) Top panels: two sub-band spectra (left) 

- Zeeman splitting in the absence of SOI, (center) -BE Bsoi , 
angle between magnetic field and ^-direction y?=0; (right) 

- arbitrary angle —7r/2<(^<7r/2. Bottom panels: magnetic 
field dependence of reflection amplitudes \r yiy \ for the spec¬ 
tra shown on top panels. For illustration we performed all 
calculations with model barrier V(y) = Vo exp(— \y\/ Lqpc), 
koFLQPc = 3.6 and hight of the barrier Vo is tuned to get 
7*0 =0.1 (see details in the text). Insert shows relative orien¬ 
tation of B and Bsoi- 


the orientation of magnetic field B with respect to y- 
axis (see Fig. 2 left bottom panel insert). 

The kinetic energy term of (|9| is given by 


n 0 


k 2 
2m* 

ctRk — i'yBe~' lcf 


apk + i'yBe 
k 2 
2m* 



( 10 ) 


The Hamiltonian Ho is promptly diagonalized in k- space. 
The eigenvalues (spectra) describe two sub-bands (v = + 
and v =8 —) split both horizontally due to the SOI and 
vertically due to Zeeman effect (Fig. 2): 


k 2 /-^^- 

E±(k) = —- ± y (a R k) + ( 7 B) - 2a R "fkBsm.^. 

( 11 ) 

(we use the short-hand notation 7 = gns/2). In the pres- 
ence of both fields there are two sub-bands while the spin 
polarization changes continuously as one moves from one 
Fermi point to the other along each sub-band. We as¬ 
sume that magnetic field is applied parallel to the plane 
of the 2DEG to decouple it from the orbital degrees of 
freedom and concentrate on the Zeeman effect only. The 
angle (p characterizes the orientation of B with respect 
to the axis of 1-d transport (y). The spectra for cp= 0 
- perpendicular orientation of B and Bsoi describe the 
situation when magnetic field B is oriented along direc¬ 
tion of the transport. If B is larger than Bsoi, the most 
important effects on thermoelectric transport are due to 
the Zeeman splitting of two sub-bands (Fig 2 upper left 
panel) 17 . In that limit the effects of B are associated 
with breaking of the channel symmetry, |r^|^|r^| (Fig 2, 
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lower left panel) which is crucial for the fate of NFlP^l 
For the case B < Bsoi we distinguish two cases: i) ip = 0 
(Fig. 2, central panel) and ii) —7r/2<(^<7r/2, ip / 0 (Fig. 
2, right panel). Since the orbital effects are negligible for 
small magnetic fields if B < Bsoi , the case i) p = 0 (Fig. 
2, central panel) can also be realized when magnetic field 
is perpendicular to the plane of 2DEG. Besides, the the¬ 
ory discussed in the paper i s also applicable when both 
Rashba and Dresselhaus SOI 21 22 are present. The trans¬ 
port coefficients for generic situation of in-plane 5-field 
are fully determined by the angle 4>=(/?o —ip between B 
and Bsoi > where <po (ip) are the angles between Bsoi 
(B) and axis of Id motion respectively. 

The eigenfunctions of Hid are momentum-dependent 
spinors ( y ) = e lk ' y Xv (k) 


X± (k) = 


_ J_ f ±ie~ i * ( - k '> 


V2\ 


( 12 ) 


where 


the first order of the backscattering potential are propor¬ 
tional to the amplitude of B (Ref. [55J 


k+/-l = n) 


kpF 

koF T $ 



cos <p, 


(17) 


where &of is the Fermi momentum at zero splitting 
((5=0), S=m*aR , ro oc \V(2koF)m*a\ < 1 is a coeffi¬ 
cient characterizing the transparency of the barrier. The 
intra-band scattering is completely suppresses for ip=7 r/2 
since the angle r d(kF±)=±i r/2 and the eigenfunctions do 
not depend on B. The intra-band reflection amplitudes 
|r+| = |r_|=0 while inter-band (\r±\,\r T \)^0. Thus, for 
ip—7 r/2 we have only two non-zero reflection amplitudes 
\r ±| and \r T \ and therefore the thermoelectric transport 
can be described by equations of Ref. |T7] if replacing 

KM r ±l and KM^I- 

The inter-band scattering amplitudes (blue dashed ar¬ 
rows on central upper panel of Fig. 2) for the case 
—7r/2<(/9<7r/2 are given by 


d (k) = arctan 


a^k — y5 sin ip 
y5 cos ip 


(13) 


The four reflection amplitudes in the first order of the 
backscattering potential are determined by 2k f momen¬ 
tum transfer and given by the matrix elements of V(y) 
in spinor basis 4/^. The diagonal matrix elements 


I v (fc£ + k^_) $(&£_) \ 

lwl_ v 0F \ 2 J 

(14) 

characterize the intra-band scattering (we shall use the 
short-hand notations ^^1 = 1^1 below). Here k^ + > 0 and 
kp_ <0 stand for the right and left Fermi points of a 
sub-band /i respectively, vq F^(m*a) _1 originates from 
the high energy cutoff, a is a lattice constant. The off- 
diagonal matrix elements (|rq | = |r±| and |r ^| = |r T |): 


| |_ V ( fc F+ ( d ( fc F+) ( fc F~) \ 

1 ^ Vo F [ 2 J 

(15) 

describe the inter-band scattering. 

The backscattering Hamiltonian 0 in the basis of 
eigenfunctions © casts the following form: 


Hbs = v F XP 1 7>l + h - c ] ■ ( 16 ) 


A fiu 


Let us analyse various limits of the back-scattering corre¬ 
sponding to different orientation of the in-plane magnetic 
field B in the regime of strong interplay with the effects 
of SOI. 

For the most generic case of interplay between SOI 
and Zeeman magnetic field there exist four independent 
scattering processes resulting in four different reffection 
amplitudes (Fig. 2, right lower panel). The reflection am¬ 
plitudes for the intra-band scattering |r+| and |r_| (black 
and red dashed arrow on Fig. 2, right upper panel) in 


\ r ±/p =r 0 


It b 


B sin ip \ 
Bsoi J 



Here coefficients (6, c(ip)) ~ 1 depend on the geometry 
of the QPC. One can see that for ip = 0, the scattering 
term linear in B (linear Zeeman effect) disappears and 
additional symmetry \r± \ = |r T | emerges (Fig. 2 central 
panel). The reflection amplitudes depend on magnetic 
field quadratically (quadratic Zeeman effect). Thus, the 
scattering Hamiltonian in that case contains three inde¬ 
pendent scattering parameters. 


V. EFFECTIVE MODEL 

We recapitulate briefly the main steps of the deriva¬ 
tion of transport coefficients (for details see Appendix A): 
i) we bosonize the 1-d Hamiltonian ( |5]|7| ) using a stan¬ 
dard approach 31 32! . The effective bosonic Hamiltonian 
gives us a boundary sine-Gordon (BSG) modef 32 ! with 
four different backscattering amplitudes. The high-T re¬ 
sults are obtained by perturbative expansion (in reflec¬ 
tion amplitudes) around the strong-coupling fixed point 
of the model, ii) The non-perturbative results in the 
low-T regime are obtained by re-fermionization proce¬ 
dure through the mapp ing t he BSG model onto the ef¬ 
fective Anderson model 12 15 17 . iii) The effects of the 
Zeeman field at the QPC resulting in TRS breaking 
caused by the asymmetry of reflection amplituded^ are 
accounted by a magnetic field dependent resonance width 
T at the CB peaks 17 . The resonance width T in the pres¬ 
ence of Zeeman field remains finite for whole range of 
the gate voltages, cuts the temperature-dependent loga¬ 
rithm and therefore restores FL properties. The width 
T (Refs. [12] [14 1 and 32]) is attributed to a single local 
Majora na mo de interacting with a single mode of chiral 
fermions 12 ^ in the theory containing only two (intra¬ 
band) scattering processes. The scattering on quantum 
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impurity in the presence of SOI and Zeeman fields in¬ 
volves four ”2 kp” processes which can be accounted for 
by two local Majorana modes interacting with four modes 
of two species of the chiral fermions. As a result, two dif¬ 
ferent resonance widths enter the transport coefficients. 
The interplay between two widths associated with inter- 
and intra-band processes leads to remarkable effects in 
thermoelectric transport. 

The effective Anderson model which describes a hy¬ 
bridization of two local Majorana fermions 771 and 772 with 
two species of c hiral fe rmions (see Appendix A) is a direct 
generalization of 12 14 16 17 for a case of interplay between 
Zeeman and SOI fields: 


r 

V) = / 

J — c 


H T (t) = I dk 

/ —OO 


Y ■ V F) C a,k C ^k 

_a=l,2 

-V 2 (w ST ( 4)771 (ci >fe -c\ k ) - iuj aT (t)rii(ci ik +c\ k ) 

+u msT (t)ri2(c2,k ~ 4,fe) “ iw mar (t)m(c2,k + 4,*: 

where following Ref. fl 6 l we define 

^s/a r(f) ^ s/a ’ fs/a r(f)y 

^ms/ma r(f) ^ ms/ma ' fs/a r(f)i 


(19) 


( 20 ) 


with 


^ s/a 


a 


ms/ma 


vfE c (P 

27r 3 

■ v F E c e c | 

27T 3 


lk+1 ± Ml, 


M ± Ml, 


( 21 ) 


and the time-dependent functions 


fs/a r = (-l)" T(t) Re/Im [exp {* (5xAt) - 7rJV)}] . (22) 


Function £y T (£) describes the deviation of the phase of 
the charge mode mean value from irn T (t) (Ref. Il 6 |) : 

&Xr{t) « |T (cot[ 7 rT(t - r)] - cot[ 7 rTf]), (23) 

here TV is a function of a gate voltage V g : N is integer 
in the Coulomb blockade valleys and half-integer in the 
Coulomb blockade peaks, n T {t) = 0(t)0(r — t) {6ft) is a 
step function) and C « 0.577 is the Euler’s constant. 

The original Matveev modeP^ corresponds to the 
case Cl a =fl rna =Cl rns =0 and describes a single Majo¬ 
rana fermion r\\ coupled to the odd combination of cre¬ 
ation/annihilation operators of chiral fermions. Contrast 
to conventional An derso n model which preserves U{ 1) 
symmetry, the model 12 is characterized by Z 2 symme¬ 
try instead. As a result, the NFL properties associated 
with the two-channel Kondo physics emerge. The NFL 
behaviour of the two-channel Kondo model are attributed 
to overscreened regime realized when the number of or¬ 
bital channels A f exceeds twice a spin of a quantum impu¬ 
rity. As it was shown in Ref. il7] Zeeman in-plane mag¬ 
netic field restores the U( 1) symmetry through appear¬ 
ance of non-zero $d a and therefore leads to the restoration 


of the FL behaviour characteristic for the single-channel 
fully screened Kondo model in both thermodynamic 15 
and transport 1 ^ coefficients. The Kondo temperature TV 
is of the order of the charging energy E c (see Appendix 
A). The effective Hamiltonian (19) describing scattering 
in the presence of both Zeeman and SOI fields has a 
structure of two copies of the two-channel Kondo model 
where coupling constants oq r depend on the magnetic 
field. Thus, when all reflection amplitudes are different, 
the model (19) is characterized by generic FL proper¬ 
ties. However, if accidental degeneracy fine-tuned by the 
orientation of in-plane magnetic field appears, one of the 
non-identical copies of the two-channel Kondo model pre¬ 
serves the NFL properties. 

The effects of interplay between the backscattering at 
the QPC and Coulomb interaction in the QD can be ac¬ 
counted by the correlator K(r ) = (T r F(r)F\ 0)) (see 
details in Refs. M and 113 . The operator F(r) accounts 
for the weak charge quantization in the mesoscopically 
Coulomb blockaded QD. Following Ref. [I6|we define the 
charge of QD Q = e(h T +h^), where h T is an integer val¬ 
ued operator which commutes with the annihilation oper¬ 
ator of the electron in the dot at the position of the source 
Ref. [HI Since by definition [F(r), n r \ = F{r ), the role of 
operator F{r) is to account for the effects of interaction in 
QD: x{t)=F{t)^q\{t) where Ta and Tqa correspond 
to interacting and non-interacting left/right fermions, re¬ 
spectively. Thus, the dressed Green’s Function (GF) 
G(t) = — (T t Ta(t)T^(0 )) and free fermionic Green’s 
function Go(r}=— (T t T 0 a(t)T^(0))=— ttuqT/ sin(7rTr) 
are connected 16 by simple relation G{r) = K(r)Go(r) 
(here v 0 is a density of states in QD). The transport 
coefficients of the model are determined by the Green’s 
function Q (see next Section). 

In order to compute the Green’s function Q (or, 
equivalently, compute the correlator K ) we define 

the operator U T =(—l) d d where d— (771 +^772 )/v/2 and 
d^ = {Vi~ ^ 772 )/\/2* We keep notations of the Matveev and 
Andreev work Ref. M for U T =2ir]2 771 - Notice, that the 
’’spin” and charge are completely disentangled in the cor¬ 
relator K(t)=K c (t)K s (t). 

While iV c (r)=7r 2 Te -C /(2F c |sin(7rTr)|) (see Ref. 14 
and m for details of calculations), the K s (r) in zero-th 
order in H r (t)—H r= o(t) is defined by the correlator 


Ki°\r) = (T t U( T )U( 0)) 0 . (24) 

Here (...)o denotes an averaging with © taken at r= 0 . 
(Notice obvious correspondence rft—>cr x /y/2, r) 2 ^cr y /x/2 
and 2 ir] 2 T]i^cr z: where cq are spin 8 = 1/2 operators 
(' i=x,y,z ). Therefore, K°(T)=(T t a z (T)a z (0)).) 

The first non-vanishing order in H r (t)—H r= o(t) cor¬ 
rection to the correlator K s (r) is given by: 

r i/T 

Ki 1 \r) = - dt(T t H' T (t)U (t)U(0))o, (25) 

Jo 

where the Hamiltonian H' r ft ) has the form: 

H' T (t ) = -2 iSx T (t) (Q s sin(7r.ZV)?7i9L - cos(7rjV)?7iCi 
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+Q ms sin(7riV)?72?2 - Yl ma cos(wN)r] 2 ( 2 ) • (26) 

Here we define four additional Major ana fermions 
(a=l,2) through a ^-Fourier transform of the even/odd 
combinations of creation/annihilation operators of two 
species of chiral fermions c a ^ taken at the position of 
the quantum impurity y= 0 : 

/ °° ^ f°° / \ 

dkCak = J y^a,k T ’ 

/ °° ^ f°° / \ 

dk<; ak = —J= J dk (c*,* - 4 )fe J . (27) 


In order to compute the correlators (24, 25) we apply 


the Wick’s theorem to the product of even number of 
fermions and express the result in terms of the products 
of the single particle GF’s. The imaginary time (Mat- 
subara) GF’s form a 6 x 6 matrix with 21 independent 
components (6 diagonal and 15 off-diagonal): 


Here we denote the resonance widths associated with 
symmetric and antisymmetric combination of reflection 
amplitudes as: 


r si ms = l/ms COS 2 (ttN)4tt/v F , 
r a/ma = ^l/ma 5 ^ (nN)4ir/ v F . 


(30) 


Two resonance Kondo widths entering the transport co¬ 
efficients are given by: 


t a r„ 


t b = t s + r a . 


(31) 


Notice that ten GF’s, namely Gj^, G^, and G^ do 
not depend on the local Majorana fermions describing 
the quantum impurity. These GF renormalise the 
correlations between the conduction electrons, but do 
not enter the Eqs. (24, 25). Another five GF’s allowed 


by the symmetry of the Hamiltonian (19) do not appear 
in the equations (24j 25) due to specific form of fermionic 
correlations in the Hamiltonian H' r (t). 


<ys( r ) = -( t t^(t)^(0)), G c Ji(t) = ~(T t Ut)U0)), 
= “(TV^(r)^(O)), G%{t) = -(T T C A1 (r)^(0)), 
GjZir) = -(T T C M (r)^(0)), G^(t) = -(T T ^(r)^(0)). 


The GF’s of quadratic Anderson-type Hamiltonian (19) 


can be found exactly (e.g. by solving equations of motion 
for the Majorana fermions). For computing the correla¬ 
tors ([24, 25) we need only 6 GF’, namely, two diagonal 


local Majorana’s GF (here R denotes the retarded GF’s): 

1 / \ 1 


GV R {e) = 


') ^ 22 ,i?( e ) - 


6 + iT B 

and four off-diagonal hybridized GF’s: 

I ( \ _ sm(7rN)27r/v F 

^n,i?( e ) = 

^22 ,r( 6 ) = 

^22 ,r( 6 ) = 


e + iT A 


(28) 



e + iT b 

Tl s COs(7tN)27t/vf 


e + iT b 

ma 

sin(7rV)27r /vp 


6 + iT A 

^ms 

COs(7tN)27t/vf 

e + iY A 


(29) 


VI. TRANSPORT COEFFICIENTS 

The thermoelectric coefficient Gt and electric conduc¬ 
tance G 


Gt = 
G = 


i7T 2 G L T 
2 e 
ttG l T 


L 


L 


sinh(7rTt) 
_oo cosh 3 (7rTt) 

1 


K 


— + ii) dt, (32) 


9 -A ( —— T it 
_oo cosh 2 ( 7 rT£) V 2T 


dt 


(33) 


are here calculated by accounting for interaction effects in 
the QD through the correlator K(r) defined in the previ¬ 
ous Section. The conductance of the barrier between the 
source and QD Gt — 27re 2 z/o^L^| 2 is expressed through 
Fermi’s golden rule as function of the density of states 
(DoS) of the source i /£, the DoS of the QD and the 
weak tunnelling amplitude \t\ (Ref. l30|) . 


The correlator Ks°\l/(2T) + it) defined by (24) is an 
even function of time. This correlator determines the 
behaviour of the electric conductance G, but does not 
contribute to Gt- 


(0) _ G L r A r B e~ c (y a y b 

G ~ Fg [t^ 


(34) 


The equation for the thermoelectric coefficient Gt is given by an odd function Ks 1 \l/(2T) + it) defined by (25): 

| r + r. 


( 1 ) _ G L r A r B sin (2ttN) 


r<\ L ) — _ 

Ajr rji - 


6e7 t 2 E c 


Y b 


In 


Eg 

T + Y b 


FI E* ^ 
4* ’ 4 


y a 


In 


Eg 

t + y a 


fi hi 

4 ’ 4 


(35) 


The ratio of Gt and G defines the thermoelectric power: 


S = 


16e c sin (2? tN) T 


2>ei r 


\r+r-1 


In 


E c \ F 


+ GM ln 


r B \t + y b J F G (^,Kf) y a \t + y a J 


Ec \ iW, 1 #) 


(36) 
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The functions Fq and F universally depend on the ratio of the resonance Kondo widths r^, T# and the temperature: 


/ OO roo 

/ 

-oo J — 


F g 


—oo J —oo 



\z + z'y 

2 +7T 2 

1 

(z') 2 + X 2 

[z 2 + y 2 } cosh (f) cosh (f) cosh (HG 


(37) 


F (x, y) 



dzdz' 


z(z + z') 

\z + z'f + 7 T 2 " 

1 

[z 2 + X 2 ] 

(z') 2 + y 2 

cosh (|) cosh (y) cosh 


(38) 



FIG. 3. (Color online) Left panel: Gate voltage and mag¬ 
netic field dependence of Y a /E c : red - Y A , blue - Yb- Right 
panel: a ”phase diagram” - four main regimes of the ther¬ 
moelectric transport are inside the green, red, blue and ma¬ 
genta domains; A - perturbative NFL, B - weak partial NFL, 
C- strong partial NFL, D - non-perturbative FL (see de¬ 
tails in Section VII). Domains boundaries are defined by the 
crossover conditions Y A (B,N)=T and Yb(B, N)=T . For all 
plots aR=0.15vF, <p=7r/4 , t*o=0.1, /cofLqpc= 3.6, T=0.05F? c . 


VII. RESULTS AND DISCUSSION 


A. Four main regimes of thermoelectric transport 


The non-perturbative equation for the resonance width 
T related to the inter-band scattering fl30|3lt demon¬ 
strates a weak dependence of Y on the magnetic field 
away from the CB peaks (see Fig.3 left panel): 

r A oc r 0 [(1 - A^) cos 2 (7 tN) + A\] , (39) 


where Yo=rQE c , (F?, cp)=b(B/Bsoi) sin tp and A\ = 
A a + 2 c((p)(B/B S oi) 2 - 

In contrast, the resonance width Y associated with the 
intra-band scattering d30|3lt strongly depends on B at 
all gate voltages: 


Y b oc r 0 [(1 



The =r B (N= |) ocr 0 • (R cos (f/B c ) 2 is a mini¬ 

mal resonance width, AB=S/koF and B C ~D, where 
D~(ra*a 2 ) -1 is the bandwidth. Thus, B c corresponds 
to the field strength that is necessary to reach full spin 
polarization of the conduction channel. 

Varying the temperature, gate voltage, amplitude and 
direction of the magnetic field one can achieve four dif¬ 
ferent regimes of thermoelectric transport (Fig. 3 right 


panel): 

A) (r A ,r B )<T - fully perturbative NFL regime. While 
Y b is gapped and the gap is Y a could be gap¬ 

less if the gate voltage is fine-tuned to the positions of CB 
peaks TV— >1 /2 and <^—»0. The TP ( 34|[36] ) demonstrates 
fingerprints of weak NFL logarithmic behaviour: 


S oc 7 q In ( ^ ) sin(27rTV). 


(41) 


B) Y b <T<B a - perturbative in Yb/T and non- 
perturbative in Y a /T (see (34 36])) . This regime can be 


reached either by fine-tuning the gate voltage away from 
the CB peaks or by tuning the direction of Zeeman field 
to be parallel to SOI in order to suppress the intra-band 
scattering: 


S oc 


\r+r -1 In ( + \ r ± r T I jL ln 


sin(27r N). 
(42) 

The weak NFL effects are manifested in the TP log- be¬ 
haviour originated from the intra-band scattering. The 
inter-band processes result in appreciable FL corrections 
to the TP ( 34p6 ). The NFL effects are weak since the 
amplitude of intra-band scattering is small at B < Bsoi • 
C) Y a <^T<^Xb - perturbative in Y A /T and non- 
perturbative in Y b /T (see 34][36] ) . This regime is 
achieved in the vicinity of CB peaks and characterized 
by strong NFL effects due to a weak magnetic field de¬ 
pendence of \r±\ and \r T \ protected by SOL Thus, by 
fine-tuning the orientation of magnetic field perpendic¬ 
ular to SOI <p=0 one can controllably protect the NFL 
behaviour of TP in the regime B<Bsoi • The magnetic 
field controlled gap associated with the intra-band scat¬ 
tering weakly depends on the gate voltage and results in 
small FL corrections to TP (compared to NFL effects): 


S oc 


\ r + r - I ln (f~) + l r±r= F I ln CY 


sin(27rTV). 
(43) 

D) T<^(r^, Yb) - FL non-perturbative regime. The NFL 
logs associated with the intra- and inter- band scattering 
processes are cut by the corresponding resonance widths 
( p4p6| ): 


S oc T 


\r+r-\ 


ln 


Ec 

r B 


k± r =Fl ln ( 

I\4 \I\4 


sin (27 tN) 
(44) 

















































FIG. 4. (Color online) Top left: eSmax as function of 
B/Bsoi for different angles p at T=0.001F C and oh=0.15?;f, 
from top to bottom (£?=7r/12,7r/6,7r/3,57r/12. Insert: eS(N) 
for (£>=571-/12 for B—D.'bBsoi (black), B—Bsoi (red), 
B=1.5Bsoi (green). Top right: eSmax as function of the an¬ 
gle p for different amplitudes of B/Bsoi at olr=0.15vf and 
T = 0.001F? C from top to bottom: B/Bsoi =0.6, 0.8,1.0,1.2.. 
Insert: eS(N) for (£>=57r/12, B=0.5Bsoi , T=0.001 (black), 
0.01 (red), 0.1 (blue). Bottom left: eSmax as function 
of B/Bsoi for different T/E c at aR=0.15vF and p—7 t/6 ; 
T/E c — 10 _1 - blue (regime A), 10 -2 - green (crossover A->C), 
10 -3 - red (B^Da-C), 10 -4 - black (BaD) . Bottom 
right: eSmax as a function of T/Vff 171 (r™ n =FB(l/2)) with 
clr=0.15vf, B=0.5Bsoi , for p=5tt/12 (black), p— 7r/3 (red), 
p—7 t/6 (green), p— tt/12 (blue), (ro=0.1 and /cqfTqpc=3.6). 


The TP is linear function of the temperature. However, 
the coefficient in front of T strongly depends on both gate 
voltage and magnetic field. 


B. Possible experimental realization and ’’smoking 
gun” predictions 

Choice of a material. We suggest to use narrow-gap 
semiconductors, e.g. InSb or In As for the observation 
of the NFL fingerprints in the quantum transport. Both 
materials are characterized by large bulk ^-factors, e.g. 
\g\ ~ 10 in InAs (see Ref. |34]) and \g\ ~ 50 in InSb 
(see Ref. |35j. The domain of parameters favourable 
for the observation of the NFL regime is defined as: 
S 0 < T < 7 B S oi < E c < e F , where S 0 ~ l/(n 0 Vqd) 
is a single-particle mean-level spacing in the QD of the 
size Lqd and ’’volume” Vqd‘ Sq D ~ h 2 /(m* Lq D ) and 
£ 3D ~ (&fTq £>) _1 * h 2 /(m*L q D ), m* is an effective 
mass of the carrier in a semiconductor, other parame¬ 
ters are defined in the previous Sections. The condition 
Lqd < lsoi allows to disregard the effects of the SOI 
in the QD (Ref. [29]), while the condition Lqpc < l m fp 
defines a ballistic regime of quantum transport through 
the QPC (here lmfp is an elastic mean-free path). 


According to Ref. [35] the parameters for InSb QD 
- QPC are as follows: hap ~ O.leV • A A 0.2eV • A, 
hoi = h/(m*ap) ~ 200nra A 400nm, \m*\ ~ 0.015m e 
( m e is electron’s mass), Esoi = m*a\/2 ~ 50/ieF, the 
typical charging energy E c ~ 1 meV and typical Fermi 
velocities are hvp ~ 0.5eV • A A leV • A , while the 
mean level spacing <So < 10/aeV for Lqd ~ Ison The 
mean-free path of the QPC of a width cIqpc ~ 10 nm is 
lmfp ~ 300 nm A 1/am. Typical parameters for the InAs 
QD - QPC are not much different!^ |m*| ~ 0.03ra e , 
hap ^ 0.05eR • A -F 0.3eF • A, hoi ~ 200 nm A 1 fim 
and l m fp ~ 300 nm A 1/im. Therefore, if we assume that 
Lqd ~ Lqpc ~ 300 nm A 500nm, our predictions could 
be verified at magnetic fields B < 500 mT and tempera¬ 
tures T ~ lOOmKA300mK for typical densities of 2DEG 
n^ deg ~ 10 ll cm~ 2 A 10 12 cm~ 2 . This estimation for pa¬ 
rameters is taken from available literature (to our best 
knowledge), but may vary due to anisotropic character 
of the ^-factor which in turn depends on external mag¬ 
netic field 34 and may also be strongly reduced in confined 
geometries of the nano-structures. 

Testing a Mott-Cutler law. The first important test of 
the interplay between effects of the SOI and Zeeman field 
is to verify the Mott-Cutler (MC) lav 33 at external in¬ 
plane magnetic field. The MC-law is a standard bench¬ 
mark for the FL propertieJ^. The MC-law says that 
the TP is proportional to a log-derivative of the electric 
conductance with respect to a position of the chemical 
potential (gate voltage V g ): 


T <91 nG 

E~cdN(V g ) 


(45) 


In the limit T (Taj'Lb) E c corresponding to the 
FL regime we get: 


<91nG 

dN(V g ) 


OC E r 


\r+r -1 \r±r T \ 


l (27riV) (46) 


while the TP is given by (44). Thus, a strong devia¬ 


tion from the MC-law in the FL regime at finite in-plane 
magnetic fields is a pre-cursor for the NFL behaviour 
discussed in the paper. Notice that break down of the 
MC-law indicates that there exists no equivalent clas¬ 
sic electric circuit consisting of the resistances connected 
in parallel or in series and therefore the effects of both 
intra- and inter-band scattering play an important role in 
the quantum transport. The violation of MC-law in the 
thermo-electric transport through a single-electron tran¬ 
sistor has been reported in Ref. [19] We are not aware of 
existence of theoretical explanation of this effect in the 
framework of the FL theory. 

Thermopower in the presence of the external B-field. 
The next step is to measure the thermopower of a 
prototype nano-device (Fig.l). The magnitude and ori¬ 
entation of the in-plane magnetic field can be controlled 
in a standard way by four magnetic coils (not shown in 
the picture). The TP maximum eS m ax{E) demonstrates 
a non-monotonic magnetic field dependence (strong 
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NFL) which is most pronounced when B is orthogonal 
to Bsoi (black curve on Fig. 4 top left and right 
panels.) This has to be contrasted to almost monotonic 
TP maximum behaviour (blue curves) characteristic for 
weak NFL - FL regimes. The non-monotonic behaviour 
of TP maximum as a function of magnetic field is a 
central prediction of our paper. The non-monotonicity 
indicates that the NFL regime of TP is protected by 
SOI contrast to FL-like behaviour demonstrating rapid 
decrease of TP when magnetic field increased^. Another 
indication of the NFL behaviour is attributed to th e ga te 
voltage dependence (Fig. 4 inserts). According tcP® it 
is characterized by strongly non-sinusoidal form (Fig. 4 
inserts). The TP maximum at zero field, NFL regime, 
scales according to 16 as eS max ^r 0 y/T/E c \n(E c /T) . The 
TP maximum in the FL regime scales as eS~T /T g q with 
T eE /E c =B/B c \n- 1 {B c /(B\rQ\ ) (Ref. H7|). The 5-field 
dependence of TP maximum measured at different 
temperatures (Fig. 4 left bottom panel) allows to 
distinguish between four main regimes A-D discussed 
in the previous subsection. This measurement can be 
used for identification of crossovers between different 
domains. The TP maximum depends linearly on T in 
the FL regime for T<T™ in =T B (N=l/2) (Ref. Ej). This 
regime holds for ip-t tt/2. In contrast to FL regime, the 
temperature dependence of TP maximum pronouncedly 
departs from the linear behaviour (see Fig. 4 right 
bottom panel) when cp is detuned from tt/2. We suggest 
to test experimentally this effect as a benchmark for the 
NFL physics. 

VIII. SUMMARY AND CONCLUSIONS 

We have demonstrated that the theory describing scat¬ 
tering of electrons characterized by two orbital degrees of 
freedom on a spin s=1/2 quantum impurity (two chan¬ 
nel Kondo model) is strongly modified in the presence 
of both appreciable spin-orbit interaction and Zeeman 
splitting. It is shown that, on the one hand, the lack 
of spin conservation due to SOI leads to the appearance 
of new (extra) scattering channels which potentially en¬ 
hance the thermoelectric transport. On the other hand, 
the Zeeman splitting produces non-zero resonance widths 
of Majorana modes describing the quantum impurity and 
thus suppresses the NFL effects. The interplay between 
these two tendencies can be controlled by fine-tuning the 
angle between Zeeman and SOI fields. Our calculations 
predict a strong dependence of the thermoelectric power 
on the angle between B and Bsoi and thus open a pos¬ 
sibility to control the scattering mechanism by changing 
between four, three or two independent scattering pro¬ 
cesses. While the cases of four- and two- weak back- 
scattering do favour the FL behaviour, the additional 
degeneracy in scattering amplitudes appearing for three- 
scattering models due to SOI effects protects the NFL 
behaviour for the range of magnetic fields B<Bsoi • We 
conclude therefore, that SOI can indeed protect the NFL 
against the destructive effects associated with breaking 
of channel symmetry. 


IX. ACKNOWLEDGEMENTS 

We are grateful to J.C. Egues, V. Fal’ko, L. Glaz- 
man, K. Kikoin, A. Komnik, S. Ludwig, C. Marcus, K. 
Matveev, L.W. Molenkamp and O. Starykh for illuminat¬ 
ing discussions. T.K.T.N acknowledges support through 
short-term visiting program of ICTP and Vietnam Na¬ 
tional Foundation for Science and Technology Develop¬ 
ment (NAFOSTED) under grant number 103.01-2014.24. 


Appendix A: Effective Hamiltonian 
Backscattering: from fermions to bosons 


The backscattering Hamiltonian mixes the left- and 
right- moving fermions: 


Hbs = Ef ^ |r 

XfjLV 


HV\ 


*L.( 0 )*a,„( 0 )+h.c. . (Al) 


The Hamilto nian s Eqs. (3-4) of the main text and the 
Hamiltonian (Al) can be bosonized 31 32 in terms of dual 


fields 4> v (y) and 9 u (y) satisfying commutation relations 
[01/(2/), = -5r<Wsgn ('V ~ y') / 2 (Refs. EHand|32|): 


*\AV) = 4=e iAfc ^ 

V27T a 


exp{i[-A0 1/ (y) + 6 v {y)]}{A2) 


where u\. v are Klein factors 31 32 introduced to ensure 
proper anticommutation relations between the right- and 
left- moving fermions. 

Using a standard procedure 3132 we introduce 
the symmetric (charge) and antisymmetric (’’spin”) 
dual variables 0 c ,s(2/) = [0+(2/)±0-(2/)]A/2 and 
0 c ,s{y)=[0+{y)±0-{y)}/V2 satisfying the commuta- 
tion relations [4> c / s (y), ^c/sW)} =-wrsgn (y - y’) /2 
(notice that we still refer to the antisymmetric in the 
band index bosonic field as ’’sp in”). We rewrite the 
backscattering Hamiltonian (Al) in terms of the charge 
and spin bosonic fields as follows: 


2D / 

H B s =-(r s cos[v / 2</> c (O)]cos[v / 20s(O)] 

+r a sin[\/2^ c (0)] sin[\/20 s (O)] 

+r ms cos[\/2<?!> c (0)] cos[\/26> s (0)] 

+r ma sin[v / 2<^c(0)] sin[\/20 s (O)]) , (A3) 


where r s = ||r+| + |r_||/2, r a = ||r+| - |r_||/2, r ms = 
lk±| + |r T ||/2, r ma = ||r±| - |r T ||/2. 


Backscattering: Majorana fermions 

As a first step we replace the charge mode by its mean 
value averaged over fast charge degrees of freedom using 
the functional integral technique developed in Ref. [16] 
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and obtain the Hamiltonian: 


/ oo 

J[d y e a (v )] 2 + [dyMv)] 2 }dv 

~\ — fw S r(i)cos[v / 20 s (O)] + u aT (t) sin[V 2 (f> s ( 0 )] 

V V]7 \ 

+w msT (t)cos[\/26» s (0)] + w mar (t)sin[v / 26» s (0)]^ , (A4) 


where we use 


the notations ( 20|21 ) of the Section V. 


As a next step we introduce the even and 

odd combinations of the ’’spin” (aka sub- 

band) bosonic fields 4> e / 0 {y)=[(j) s (y)±(j) s (-y)]/2 , 
0e/o{y)=[0s{y)±Qs{-y)\/'/ 2. As a result, we obtain new 
chiral fields ^ 1/2 {y) = ^o/e{y)~4 > e/o{y) satisfying the com¬ 
mutation relations: [& a (y)^ a / (y r )]=^S aoi fSgD.(y — y') 
where a, a'=1,2. We define new fermionic fields 

{y) = { 1 7a/v / 27ra) exp (— i& a (y)) with a help of 
two local Majorana fermions rji=(d+d^) / y/2 and 


V 2 =(d—cft)/(iy/2) representing the quantum impurit}^. 

Finally, we integrate out the fluctuations of the spin de¬ 
gree of freedom with the frequencies exceeding E c (Ref. 
US]). This procedure is equivalent to the poor man’s scal¬ 
ing approach originally used for the Kondo problemP^ 
and leads to replacement of the bandwidth D by the new 
bandwidth Tx ~ E c . As a result, we derive the effec¬ 
tive Anderson model which describes a hybridization of 
two local Majorana fermions iji and 772 with two species 
of conduction electrons. The effective Hamiltonian (19) 


has a structure of two copies of the two-channel Kondo 
model where coupling constants c^ r depend on both Zee- 
man and SOI fields: 


iT 


r 

V) = / 

J — ( 


dk 


E ‘ V F)cl }k Ca,k (A5) 

Loi=1,2 

-V2 (w aT (4)771 (ci,fe - c\ k ) - iu a r(t)v l( c l,fc + c l,fe) 
+u mST (t)r} 2 {c 2 ,k - c\ k) - iu) mar (t)V2{C2,k + 4 ,fe) 
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